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ABSTRACT 


A  method  for  approximate  treatment  of  coll Inear  H-atom  or  proton  transfer 
reactions  AH  +  B  -*•  A  +  HB  which  involve  two  nearly  degenerate  vibrational 
states  Is  developed.  The  method  is  based  on  constructing  a  diabatlc  represen¬ 
tation  of  the  problem,  i.e.  one  In  which  the  reactant  and  the  product  vibrational 
states  are  distinguished  as  such  throughout  the  collision.  The  diabatic 
representation  is  constructed  with  the  aid  of  the  adiabatic  one  discussed  in 
an  earlier  paper. 

The  treatment  can  also  be  applied  to  symmetric  exchange  reactions,  in 
which  case  it  yields  good  agreement  with  other  results. 


I.  INTRODUCTION 


In  an  earlier  paper\  henceforth  referred  to  as  I,  a  simplified  approximate 
treatment  of  symmetric  coll  Inear  reactions  In  which  an  H-atom  or  a  proton  Is 
transferred  between  two  heavier  particles  was  developed.  The  treatment  was 
based  on  the  adiabatic  near  separability  of  the  faster  H-atom  motion  from  the 
slower  relative  motion  of  the  heavy  particles.  It  was  useful  to  employ  polar 
coordinates:  The  radial  coordinate  represents  the  translational  motion  of 
the  heavier  particles  and  the  angular  coordinate  corresponds  to  the  vibrational 
motion  of  the  H-atom. 

Due  to  the  exact  degeneracy  of  the  vibrational  motion  of  the  reactants 
and  the  products  the  adiabatic  vibrational  states  are  symmetric  and  antisymmetric 
linear  combinations  of  the  reactants'  and  products'  vibrational  states,  which 
are  exactly  decoupled  by  symmetry.  As  a  result,  the  reactive  transition  probabil¬ 
ities  between  states  with  the  same  vibrational  quantum  number  can  be  evaluated 
approximately  within  the  adiabatic  approximation. 

There  have  been  some  accurate  numerical  calculations  on  col linear  H-atom 
transfer  reactions  since  the  publication  of  I  which  provide  an  opportunity  to 
test  the  validity  of  the  approximations  used  there.  The  comparisons  of  the 
results  of  the  accurate  numerical  calculations  with  those  of  I  are  presented 
In  Sec.  2.  A  somewhat  different  approach  to  the  symmetric  H-atom  exchange, 
which  has  the  advantage  of  being  more  readily  extendable  to  the  nonsymmetric 
case  is  also  developed  and  tested  in  Sec.  3. 

The  transfer  of  an  H-atom  in  an  asymmetric  reaction  is  considerably  more 
difficult  to  treat.  Unlike  in  the  symmetric  case,  the  reactant  and  product 
states  are  no  longer  exactly  degenerate  and  the  reaction  can  occur  only  as  a 
result  of  a  breakdown  of  the  adiabatic  approximation.  Consequently,  two  or  more 
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coupled  vibrational  states  are  generally  Involved  In  the  dynamics  and  approximate 
solutions  to  a  set  of  coupled  equations  must  be  developed. 

Reasonably  good  analytical  approximations  tothe  reaction  probability  over 
a  wide  energy  range  can  be  developed  only  In  the  case  when  the  problem  can  be 
reduced  to  the  Interaction  of  only  two  coupled  vibrational  states,  the  Initial 
reactant  state  and  the  final  product  state,  namely  when  they  are  nearly  degenerate. 
This  Is  the  case  which  Is  treated  in  the  remainder  of  the  paper. 

The  procedure  used  in  the  paper  to  treat  the  above  problem  Is  the  following: 

1.  The  diagonal  elements  of  the  equations  In  the  adiabatic  representation 
are  first  obtained  using  the  procedure  described  In  I. 

2.  The  diabatlc  equations  are  then  obtained  In  an  approximate  way  utilizing 
the  above-mentioned  adiabatic  ones  (a)  for  the  case  of  strong  and  moderate 
coupling  and  (b)  (semiclassically)  for  the  case  of  very  weak  coupling  (Sec.  4). 

3.  The  transition  probability  Is  finally  obtained  from  the  coupled  diabatic 
equations  using  an  exponential  DWBA  approximation. 

The  reason  for  step  1  is  that  the  coupled  adiabatic  equations  In  terms  of 

1-3 

polar  coordinates  constitute  a  very  compact  representation  of  the 
system.  The  reason  for  going  through  step  2  is  that  a  simple  analytic  approx¬ 
imation  to  the  solution  of  two  coupled  equations  can  most  readily  be  obtained 
in  the  diabatic  representation. 

A  perturbative  (DWBA)  method  for  solving  the  diabatic  equations  in  the  weak 
coupling  limit  and  a  procedure  for  extending  the  validity  of  the  results  to 
cases  of  higher  coupling  are  then  used  (Sec.  5)  to  obtain  the  transition 
probabilities. 

Numerical  results  are  presented  In  Sec.  6. 
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II.  SYMMETRIC  H-ATOM  EXCHANGE.  COMPARISON  WITH  ACCURATE  NUMERICAL  RESULTS. 

It  was  shown  In  I  that  the  Schrodlnger  equation  for  a  col  linear  H-atom 
transfer  collision  can  be  transformed  Into  a  set  of  coupled  ordinary  differential 
equations  In  the  adiabatic  representation: 

\-I  E  JQ1J>  *j(p)  *  0  <21> 

The  symbols  are  defined  In  I.  It  was  also  shown  that  for  a  symmetric  system 
the  reactive  transition  probability  without  a  change  in  vibrational  quantum  number 

D 

P„„  can  be  obtained  from  the  elastic  collision  phase  shifts  for  the  two  lowest 
oo  r 

states  of  Eq.  (2.1).  Those  states  satisfy 


*•:«»]♦>>  -  *»>> 

[-$£■  *  *>>]*>>  ■  e»o<c) 


(2.2a) 


(2.2b) 


which  are  exactly  decoupled  by  symmetry. 


Is  given  by: 


slnMc! 


(2.3) 


where  C*  and  are  the  phase  shifts  for  eq.  (2.2a)  and  (2.2b),  respectively. 


There  have  been  some  accurate  numerical  calculations  on  col linear  H-atom 
transfer  reactions  since  the  publication  of  I,  namely  for  the  IH  +  I  and  I  +  yl 
systems^’4  which  enables  us  to  test  the  accuracy  of  the  approximation  developed 
In  I.  Fig.  (1)  shows  the  comparison  of  the  results  of  Eq.  (2.3)  and  of  the 
accurate  numerical  calculation  for  the  IHI  system.  It  can  be  seen  that  if  the 
mass  of  the  end  atoms  is  sufficiently  greater  than  the  mass  of  the  center  atom, 
the  results  are  pratically  identical  with  the  results  of  the  accurate  numerical 
calculations. 

3  R 

The  only  difference  in  the  case  of  the  IHI  system  is  that  the  PQ0  vs. 
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energy  curve  obtained  form  eq.  (2.3)  is  shifted  by  0.5  •  10  kcal  toward  the 
lower  energies  in  respect  to  the  accurate  numerical  results  .  Such  shifts  were 
also  discussed  in  I  and  are  roughly  propositional  to  the  ratio  of  the  masses  of 
the  end  and  the  middle  atom.  For  the  Iyl  system  the  corresponding  shift  can 
be  expected  to  be  ibout  an  order  of  magnitude  smaller,  which  is  most  likely  well 
within  the  limits  of  accuracy6  of  the  numerical  calculation4. 

Eq.  (2.2)  is,  as  confirmed  by  the  agreement  in  Fig.  1.,  a  very  convenient 
one  for  evaluating  the  ground  state  reaction  transition  probabilities  for  a 
symmetric  H-atom  exchange  reaction. 

This  treatment,  however,  can  not  be  directly  generalized  to  an  asymmetric 
system  which  lacks  the  symmetry  basis  for  the  decoupling.  To  facilitate  the 
extension  to  the  asymmetric  case,  another  two  state  approximation,  one  in  the 
so-called  dlabatic  representation,  will  be  developed  and  tested  in  the  remainder 
of  this  section.  Since  in  the  two  state  approximation  (Eq.  2.2)  the  set  of 
coupled  adiabatic  equations  (Eq.  2.1)  for  a  symmetric  system  is  exactly  solvable. 
Its  solution  Is  used  to  check  the  accuracy  of  the  alternate  dlabatic  approach, 
for  which  it  constitutes  the  starting  point. 
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III.  DIABATIC  TREATMENT  OF  THE  SYMMETRIC  CASE 

The  equations  in  the  diabatic  representation  are  the  equations  for  the 
p-dependent  coefficient  in  an  expansion  of  the  wave  function  in  the  form 

ij>(s;p)  =  I  iMp)  tu  (s;p)  (2.4) 

i  1  1 

where  n^(s;p)  are  s -dependent  vibrational  wave  functions7  with  a  weak 
parametric  dependence  on  p.  They  can  be  written  in  a  matrix  form  as: 


\  ~  +  1  (p)J  i  (p)  =  e  i 


(p) 


(3.2) 


where  V  (p)  is  a  Hermitian  matrix  whose  elements  are  all  scalars, 
i.e.  do  not  contain  any  derivatives, and  £(p)  is  a  column  vector  with 
elements  ^(p). 

A  set  of  coupled  equations  in  the  adiabatic  representation  can,  in  general, 
be  transformed  into  the  diabatic  representation  by  finding  a  transformation 
matrix  T  acting  on  the  solution  wave  function  vector  ♦: 

*  T  &  (3.3) 

T  is  given  by  the  solution  of  the  differential  equation: 

p  T  «  P  T  (3.4) 

where  p  is  the  momentum  operator  p  ■  -1h^  and  P  is  the  kinetic  coupling 
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matrlx  with  elements  P. .  In  the  set  of  coupled  equations  in  the  adiabatic 

•  J 

representation.  The  appropriate  boundary  condition  for  the  equation  in  the 
diabatic  representation  is  that  reduces  to  the  properly  separable  asymptotic 
translational  wave  function  at  large  p. 

In  the  particular  case  of  (2.2),  the  coupling  matrix  JP  in  (2.2)  vanishes 
because  of  symmetry.  The  solution  of  (3.4),  namely  T  *  const.,  is  then  trivial, 
the  value  of  the  constant  matrix  T  being 


T  = 


1 

~n 


i  r 
.i  -i . 


(3.5) 


since  at  large  p  <f>s  and  <f>a  are  symmetric  and  antisymmetric  linear  combinations 
of  the  proper  asymptotic  states.  The  transformation  (3.3),  with  T  given  by 
(3.5),  transforms  (2.2)  into  the  form  (3.2),  namely: 


.  1  _il  +  Vn  (p)  -  E 
t  2  dp*  11 


^(p)  =  -v12  *2(p) 


+v“  (0)-E 


^2(p)  =  _V2i  j(p) 


where 


V11  =  V22  =  e  ; 


c+  =  E1  +  e2 


V12  *  V21  "  e  ; 


-  £a  -  C, 

c  *  2  i 


(3.6) 

(3.7a) 

(3.7b) 


Solving  the  coupled  equations  (3.6)  numerically  would,  of  course,  give 

the  same  results  already  given  by  Eq.  (2.3).  Eq.  (3.6)  can  also  be  solved 

q 

perturbationally  yielding  what  Is  usually  known  as  the  DWBA  solution  for  the 

transition  matrix  element: 

00 

T12  =  *1  f»)  v12  <“>  »2  (0)  “0  °'8) 

0 

where  ij/?(p)  are  the  solutions  of  the  analogs  Eq.  (3.6)  with  the  r.h.s.,  i.e. 
the  coupling  V12  set  to  zero.  Linearizing  the  V^.  and  the  exponent  of  V12, 
a,  (*  -lnV^2)  around  the  turning  point,  one  can  evaluate  the  Integral  (3.8) 
analytically  to  yield  (c.f.  Sec.  V): 

o  — 

W  ■  ^2  ■  ^5  (3-9> 

Eq.  (3.9)  differs  from  the  semiclassical  result  (5.11)  in  I  only  by  the  ex¬ 
ponential  factor  in  (3.9), which  Is  in  most  cases  very  close  to  unity.  A  more 
accurate, albeit  more  complicated, analytic  approximation  to  (3.8)  than  (3.9) 
also  exists10.  For  the  present  purposes,  however,  the  simpler  expression  (3.9) 

is  sufficiently  accurate. 

The  perturbative  DWBA  expression  (Eq.  3.8)  is  valid  only  in  the  weak 
coupling  case  (i|»  “  ^.|T12|*  «  1).  It  can,  however,  be  readily  extended 
to  the  case  of  higher  coupling  (i.e.,  higher  transition  probabilities)  by 
introducing  in  an  ad  hoc  way  the  sinusoidal  dependence  of  the  transition 
probability 
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Pg0  »  s1n2(2itT°2)  .  (3.10) 

known  as  the  exponential  DWBA  approximation.^ 

The  transition  probability  vs.  energy  obtained  from  Eq.  (3.9)  and  (3.10) 

is  shown  in  Fig.  2,  together  with  the  exact  solution  of  E.  (2.2)  given  by  Eq. 

(2.3).  The  logarithm  of  the  transition  probability  vs.  energy  in  the  tunneling 

region,  obtained  from  Eq.  (3.10)  and  from  (2.3),  is  shown  in  Fig.  3. 

The  surface  used  is  a  LEPS  surface  with  the  asymptotic  parameters  being  those 

1 2 

of  the  IH  +  I  I  +  HI  reaction.  Its  potential  energy  contour  plot  is  shown 
in  Fig.  4, 

It  is  seen  from  Figs.  2  and  3  that  the  exponential  DWBA  approximation  to 
the  result  of  the  coupled  equations  in  the  diabatic  representation  yields  for 
the  symmetric  case  a  result  in  satisfactory  agreement  with  the  exact  result 
(Eq.  2.3).  The  exponential  DWBA  approximation  to  the  solution  of  a  pair  of 
coupled  equations  can,  therefore,  be  assumed  to  give  a  good  approximation  to  their 
solution. 

In  the  next  section,  a  procedure  for  obtaining  a  pair  of  coupled  equations 
in  the  diabatic  representation  for  a  nearly  degenerate  asymmetric  H-atom  exchange 
problem  is  developed. 


IV.  THE  NEARLY  DEGENERATE  H-ATOM  TRANSFER.  GENERATING  THE  DIABATIC 
REPRESENTATION. 

(a)  Moderate  to  strong  coupling. 

We  consider  the  case  where  two  of  the  adiabatic  diagonal  terms  in  Eq.  (2.1), 
one  of  which  asymptotically  corresponds  to  a  reactant  state  and  the  other  to  a 


product  state,  are  nearly  degenerate  at  large  values  of  p.  The  nature  of  the 
interaction  between  such  two  states  is  closely  related  to  that  between  two  states 
which  have  the  same  quantum  number  in  a  symmetric  system.  For  the  nearly  degen¬ 
erate  case,  the  following  closely  related  approximate  treatment  is  developed. 

First,  the  adiabatic  eigenvalues  are  calculated,  either  quantum 
mechanically  as  described  in  I,  or  if  one  prefers  less  computation  at  a 

slight  loss  of  accuracy,  from  a  semiclassical  treatment  of  a  double  minimum 
13 

potential . 

The  most  important  simplification  arises  from  the  fact  that  the 
two  nearly  degenerate  states  interact  much  more  strongly  with  each  other  than 
with  the  other  states  of  the  system,  and  the  interaction  with  the  other  states 
can  be  neglected.  The  problem  can  hence  be  reduced  to  interaction  of  two  states, 
a  case  particularly  amenable  to  an  approximate  solution. 

Upon  neglecting  all  but  the  two  nearly  degenerate  states,  Eq.  (2.1)  reduces 
to 

[-1  fr  *  ^  -E  *<!„] 

(4.1) 

["  *  t7  *  *2  '  ip7  '  E  +  Qz2]  *2<pl  =  1P2'  ^  *,<P> 

An  example  for  which  such  near-degeneracy  occurs  is  the  Cl H  +  Br  -*■  Cl  + 

lO 

HBr  reaction.  A  contour  plot  for  a  model  LEPS  surface  with  the  asymptotic 
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parameters  such  as  to  fit  the  above  reactions  is  shown  on  Fig.  6.  The  diagonal 
elements  of  Eq.  (4.1)  for  this  surface  are  shown  In  Fig.  7  with  the  dashed  lines. 

The  usual  method  for  transforming  from  an  adiabatic  to  a  dlabatlc  repre¬ 
sentation  by  diagonalizing  the  kinetic  energy  matrix  can  generally  not  be  used 
when  there  Is  a  transfer  of  a  particle  involved,  except  In  the  somewhat  trivial 
case  of  Eq.  (2.2)  where  the  kinetic  coupling  matrix  vanishes.  The  reason  is 
that  a  part  of  the  parametric  dependence  of  the  basis  set  x(s;p)»and  consequently, 
of  the  kinetic  coupling,  is  an  artifact  of  the  way  In  which  the  x(s;p)  are 
defined  and  does  not  correspond  to  an  actual  interaction  of  the  underlying 
diabatic  states.  This  problem  appears  frequently  in  electron  transfer 
reactions  and  is  known  as  the  translation  factor  problem. Another  way  of 
determining  the  diabatic  matrix  elements  V^,  V ^  and  V^2  which  does  not  make 
use  of  the  nonadiabatic  coupling  element  needs  to  be  devised. 

The  di agonal izati on  condition  for  a  two  state  diabatic  potential  energy 
matrix  can  be  written  as 


1.2 


CM 

CM 

> 

+ 

[V22  -  Vlll 

2  g 

4.  \IC 

2  -  \ 

z  j 

+  v12 

(4.2) 


or,  in  a  form  more  convenient  for  our  purposes  in  terms  of  and  V-»  as 


+ 

e— 


En  -  e22  . 


V±  s  V11  1  v22 

w  o 


(4.3) 


One  finds  from  Eq.  (4.2)  that 


V+  -  e+ 


(4.4) 


-n- 


and 


H12  “  (e’)2  ‘  (V’)2  *  (4.5) 

That  is.  the  sum  of  the  two  diabatic  diagonal  terms  equals  the  sum  of 
the  adiabatic  ones  in  the  two  state  approximation.  However,  in  order  to 
obtain  the  diabatic  interaction  matrix  element  (Eq .  4.5),  as  well  as  the 

individual  diagonal  diabatic  elements  (Eq.  4.3),  the  difference  between  the 
diagonal  diabatic  matrix  element  for  the  two  states  involves  must  be  known. 

It  can  be  assumed  on  physical  grounds  that  the  diabatic  vibrational  states, 
which  in  general  can  be  identified  with  a  particular  configuration  throughout 
the  collision,  are  associated  with  one  of  the  wells  of  the  double  well  potential 
(Fig-  6).  Eigenstates  localized  in  a  single  well  can  be  defined  in  a  multitude 
of  ways:  These  eigenvalues  will  generally  be  somewhat  sensitive  to  the  choice 
and  cannot  be  readily  identified  with  the  diabatic  states' eigenvalues.  The 
difference  between  the  single  well  eigenvalues  is  substantially  less  sensitive 
to  that  choice  and  will  be  assumed  to  approximate  fairly  closely  the  difference 
between  the  diabatic  eigenvalues.  Accordingly,  we  make  the  approximation 


where  E°  and  E^  are  suitably  defined  single-well  eigenvalues.  Both  the 
diagonal  and  the  off-diagonal  elements  of  the  diabatic  potential  energy  matrix 
can  then  be  evaluated  using  Eqs.  (4.3)  through  (4.6). 
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The  single  well  eigenvalues  E°  and  E2  in  (3.7)  ,  which  now 
define  Vjj,  Vgg  an<1  vi2  vla  ^qS-  (^-3)  throu9h  (4.6)  are  strictly  defined 
only  if  the  E°  lie  below  the  barrier  top,  which  seems  to  limit  the  present 
treatment  to  tunneling  energies.  However,  the  V”  in  (3.7)  is  a  very  weak  function 
of  p,  which  makes  it  possible  to  simply  extrapolate  V”  and  thus  determine  Vjj  and 
V 22  at  smaller  values  of  p,  where  E°  and  E°  are  difficult  to  determine  directly. 
This  extrapolation  enables  one  to  extend  the  treatment  to  energies  throughout  one 
full  cycle  of  increasing  and  decreasing  P. .  in  the  transition  probability  P..  vs.' 

i  j  *  j 

energy  curve. 

The  resulting  Vjj  and  V12  versus  p  for  surface  IV  is  shown  in  Fig.  7  with 
the  dashed  line.  Log  VJ2  vs.  p  shown  in  Fig.  8.  The  right  end  of  the  curve  for 
log  V12  denotes  the  limit  of  the  reliability  of  the  small  numbers  V12  obtained 
as  a  difference  between  two  much  larger  numbers,  (Eq.  4.5).  A  much  more  reliable 
treatment  for  this  region  of  p's  (large  p's), based  on  the  semiclassical  treat¬ 
ment  of  a  double  minimum  potential, is  given  next. 

(b)  Weak  coupling,  semicalssical  evaluation  of  Hj2. 

The  interaction  matrix  element  V^2  is  evaluated  in  (4.5)  from  the  difference 
between  (W")J  and  (V’)1,  which  in  turn  are  given  by  differences  between  the  two 
adiabatic  or  diabatic  levels.  In  cases  when  the  barrier  between  the  two  wells  is 
very  large,  i.e.,  Vj2  is  very  small,  this  evaluation  of  a  small  number  as  a 
difference  of  two  large  ones  can,  and  does  lead  to  a  serious  loss  of  accuracy. 
Another  more  direct  approach  to  evaluating  Vj2  which  circumvents  the  above 
problem  in  the  lowV12  limit  is  next  obtained  from  the  semiclassical  treatment 
of  the  eigenvalues  of  a  double  minimum  potential. 

The  eigenvalues  of  a  double  minimum  potential  which  lie  below  the  barrier 
top  can  be  obtained  semicalssically  by  finding  the  energies  Ej  or  E2  at  which 
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the  semlclassical  double  well  eigenvalue  condition  (4.7) 


Is  satisfied.13*16,17 


tan 


(E)  -  f  -  ^3  tan  (>2(E)-f  -  £  J  -  F2 


(4.7) 


The  symbols  In  Eq.  (4.7)  are  given  by 


(4.8) 


(4.9) 


(4.10) 


where 

Pi  “  *  {2[E1  -  V(s;p)3}Ifi  (4.11) 

and  the  turning  points  sj  In  Eqs.  4.9  and  4.10  are  marked  on  the  cut  through  the 
surface  at  p  *  const.,  shown  In  Fig.  8. 


? 
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On  the  other  hand.  If  the  wells  are  considered  separately,  one  can  define 
distorted  single  well  eigenvalues  Influenced  by  the  presence  of  the  barrier,  but 
for  which  r(E)1nEq.  (4.7)  is  assumed,  purely  for  this  definition  to  be  zero,  l.e. , 


2n.  +  1  . 

a1<E1>  *  — 2 - "  +  I 


(4.12) 


where  n^  Is  the  vibrational  quantum  number. 

Since  one  expects  that  each  of  the  two  E^  would  be  close  to  the  E?,  one 
can  simplify  Eq.  (4.7)  by  expanding  it  around  one  of  these  E^'s.  One  obtains 


tita  -  E°)  tan  [ct2(£)  -  J  -  £}  *  r2  (E,) 
at 


(4.13) 


Keeping  in  mind  that  the  E^  are  nearly  degenerate,  one  can  further  simplify 
(4.13)  by  expanding  the  tangent  about  E°  to  get  for  the  shift  due  to  the  presence 
of  the  other  well.  Thereby,  we  have 


ft  Tiv.  flv«  n 

E,  -  E°  =  r7  e2<e,> 

tz  -  t1 


(4.14) 


where  is  the  local  frequency,  or  hv.  reciprocal  density  of  states,  given 


,  .  i  < 
1  Hrvj" 


(4.15) 


Similarly, 


e2  -  •  ^1  *v2  r2  (e2) 


(4.16) 


E  -E° 

tl  t2 
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Subtracting  the  two  level  shifts  for  the  eigenvalues  for  the  left  and 
right  well,  and  approximating  E.  -  E°  by  AE  +  AE°  one  gets: 


where 


AE  -  AE° 


flV1  flVg 
AE  +AE° 


(4.17) 


(E2  =  ^1 ) * 


AEW 


=  —  (E° 
2  1  2 


e?),  r1  =  r(E1).  r2  =  r(E2)  (4-18> 


Identifying  the  difference  of  the  two  double  well  eigenvalues  AE  in 
Eq. (4.5) with  the  difference  of  the  adiabatic  eigenvalues  e“  and  the  difference 
of  the  distorted  single  well  eigenvalues  aE°  with  that  of  the  diabatic  eigenvalues 
V~  and  rewriting  Eq,  (4.5)  as 

AE  -  AE°  =  - ^-r  (4.19) 

AE  +  AE° 


and,  comparing  (4.16)  and  4.17),  one  obtains 


fiv2 


2 

r  1 


(4.20) 


which  In  the  limit  of  a  symmetric  system  (vj  ■  v2  -  vJt  Tj  «  r2  *  r)  becomes  the 
usual  degenerate  level  splitting  result. 


V 


12 


■  for 


(4.21) 
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The  principal  condition  for  the  validity  of  Eq.  (4.18)  Is  the  near  degeneracy 
of  the  E.j,  i.e.,  that  the  spacing  of  E2  or  from  any  other  leve  in  the  system 
is  much  larger  than  E2  -  Ej,  which  is  the  basis  of  the  two  state  approximation. 


V.  EVALUATION  OF  THE  TRANSITION  PROBABILITIES 


Once  an  approximation  to  the  diabatic  matrix  elements  V^.  has  been  obtained, 
one  can,  in  principle,  solve  for  the  uncoupled  analogs  of  Eq.  (3.6) 

-  7r7^(p)  '  CE  -  v„(p)3  (5-1) 

in  order  to  obtain  the  transition  probability  via  Eq.  (3.8). 

If  v-ji(p)  in  (5.1)  has  no  multiple  turning  points,  a  good  approximation 
to  the  solution  for  the  purpose  of  using  it  in  Eq.  (3.8)  can  be  obtained  by 
solving  the  equivalent  of  Eq.  (5.1)  for  a  potential  linearized  around  the 
turning  point 

E  -  V^p)  *  F  (P  -  Pi)  (5.2) 


where  F 


dV 


1i 


P“Pj  (1*1,2),  and  p^  are  the  classical  turning  points. 


Writing  as  V°2  e”01^  and  expanding  o(p)  around  the  outer  turning  point, 
one  obtains  the  Integral  (3.8)  In  a  form  which  can  be  evaluated  analytically: 

2V12  r  «'(p-p1)  do 

T*,--^  j  A1C-B(p-P1)3  A1  C-B(p-p2)]  e  1  dp 


where  6  *  (2F)*^  and  a'  * 


(5.3) 


Using  the  Integral  form  for  the  Airy  functions,  Eq.  (5.3)  transforms  into 

-o  V12  f dx  /" dy  1/3  (xS+  r  dpe-1CB(x+y)-ia3p 

T12  *  e 


(5.4) 


where  A  is  P2~py  After  some  algebra,  one  gets 


a'  As83 


o  _  v12  Q  T2F  ‘  4a 
12  ‘  2iraF  e 


(5.5) 


where  Vj^  is  the  value  of  V^Jp)  halfway  between  the  turning  points. 


VI.  NUMERICAL  RESULTS 


The  calculations  for  the  nearly  degenerate  H-atom  transfer  were  made  on  a 

model  LEPS  potential  energy  surface  which  has  the  assn\ypotic  parameters  for  the 

Cl  +  HBr-»- C1H  +  Br  reaction.  The  Sato  parameters  were  set  to:  0.02,  0.02  and  0.0 
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In  order  to  obtain  a  surface  with  a  barrier  height  of  *21  kcal  mole 

The  adiabatic  eigenvalues  used  to  obtain  V+  via  Eq.  (4.4)  were  calculated 

19 

using  the  finite  element  method  used  in  I.  The  one-well  eigenvalues  used  to 
obtain  V"  Eq.  (4.6)  were  obtained  a)  from  potentials  such  that  the  potential  was 
set  to  a  constant  equal  to  the  barrier  top  energy  from  the  barrier  top  into  the 
other  well,  using  the  same  method,  and  b)  semiclassically  using  Eq.  (4.12).  The 
V  was  calculated  only  for  values  of  p  for  which  the  results  of  a)  and  b)  were 
in  agreement. 

Vjg  was  calculated  from  Eq.  (4.5)  and  Eq.  (4.19),  respectively.  In  the  re¬ 
gions  of  p  noted  In  Fig.  8. 

For  surface  II  the  nearly  degenerate  transition  Is  a  (0-*-2)  one.  The  tran¬ 
sition  probability  P^  was  evaluated  using  Eq.  (5.1)  and  is  given  on  Fig.  9  as 

D 

a  function  of  the  energy.  The  log  of  the  transition  probability  P02  vs.  energy 
in  the  tunneling  region  is  given  in  Fig.  10. 
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FIGURE  CAPTIONS 


Figure  1. 

Ground  vibrational  state  transition  probability  for  the  IH  +  I  -*■ 

I  +  HI  reaction.  The  solid  line  is  the  result  of  Eq.  (2.3)  and 

the  circles  are  the  accurate  numerical  results  of  ref.  3. 

Figure  2. 

12 

Ground  vibrational  state  transition  probability  for  the  model 

surface  I.  The  solid  line  is  the  exact  solution  of  Eq.  (2.2) 

given  by  Eq.  (2.3)  and  the  circles  are  the  result  of  the  present 

exponential  DWBA  approximation  (Eq.  3.10).  The  dotted  line  is  the 

DWBA  without  the  exponential  approximation. 

Figure  3. 

The  logarithm  of  the  transition  probability  vs.  energy  for 

surface  I.  The  solid  line  is  the  result  of  Eq.  (2.3)  and  the 

circles  are  the  present  result  (Eq  3.10). 

Figure  4. 
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Contour  plot  for  the  surface  I  . 

Figure  5. 

Contour  plot  of  the  surface  II19.  The  contour  energies  are  in 
kcal  mole"^ . 

Fiaure  6. 

A  cut  through  model  surface  II  at  p  ■  const.  The  solid  horizontal 

lines  show  the  adiabatic  levels  and  the  dashed  ones  the  diabatlc 
levels.  The  s]  In  the  figure  are  the  Integration  limits  for 

Eqs.  (4.9)  and  (4.10)  for  E  ■  Ej. 

Figure  7. 

Diagonal  adiabatic  (dashed  line)  and  diabatlc  (full  line)  matrix 

elements  for  model  II.  (Cl  +  HBr  -►  C1H  +  Br) 

The  log  at  vs.  p  for  model  surface  II.  The  dashed  part  of 
the  curve  corresponds  to  the  region  where  Eq.  (4.5)  becomes 
unreliable.  (Cl  +  HBr  -*•  Cl H  +  Br) 

The  transition  probability  vs.  energy  (Eq.  5.5)  for  the  model 
surface  II.  (Cl  +  HBr  -*•  C1H  +  Br) 

The  log  of  the  transition  probability  for  surface  II  in  the 
tunneling  region. 
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